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Abstract 

We propose exact vacuum expectation values of local fields for a quantum group 
restriction of the 6*2^^ affine Toda theory which corresponds to two coupled minimal 
■ models. The central charge of the unperturbed models ranges from c = 1 to c = 2, 
where the perturbed models correspond to two magnetically coupled Ising models and 
^ ! Heisenberg spin ladders, respectively. As an application, in the massive phase we deduce 



in 
o 
o 



the leading term of the asymptotics of the two-point correlation functions. 
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• • '■ 1 Introduction 

> ■ 

• 1—1 , 

>< ■ The vacuum expectation values (VEV)s of local fields play an important role in quan- 

tum field theory (QFT) and in statistical mechanics @]. In statistical mechanics the 
VEVs determine the "generalized susceptibilities", i.e. the linear response of a system 
to external fields. Furthermore, the VEVs provide all the information about correlation 
functions in QFT defined as a perturbed conformal field theory (CFT) that is not ac- 
cessible through a direct calculation in conformal perturbation theory [Q. A few years 
ago, some important progress was made in the calculation of such quantities in integrable 
(1-1-1) QFT. In ref. 0], an explicit expression for the VEVs of the exponential field in 
the sinh-Gordon and sine-Gordon models was proposed. In ref. |^ it was shown that 
this result can be obtained using the "refiection amplitude" of the Liouville field the- 
ory. This method was also applied in the so-called BuUough-Dodd model with real and 
imaginary coupling. It is known that c < 1 minimal CFT perturbed by the operators 
$12 and $13 can be obtained by a quantum group (QG) restriction of the sine-Gordon 
and imaginary Bullough-Dodd model with special values of the coupling. The VEVs 
of primary fields were then calculated. The same method was applied later to integrable 
QFTs involving more than one field. For instance, the VEV for a two-parameter family of 
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integrable QFTs ^ gave rise to the VEV of local operators in parafermionic sine-Gordon 
models and in integrable perturbed SU(2) coset CFT [jlO[. The VEV of simply-laced 



affine Toda field theories (ATFT)s is known for a long time [|TT] and the case of non 



simply laced dual pairs was recently studied in |]T2|, |I3] for which a general expression for 
the VEVs was derived. 

Such perturbed CFTs have recently attracted much attention in condensed matter 
physics, such as in the context of point contacts in the fractional quantum Hall effect 
and impurities in quantum wires |]T1|, |T3|. In such cases the property of integrability has 
provide a non-perturbative answer for experimentally important strongly interacting solid 
state physics problems [jlS, 16|. Particulary, on-shell results were obtained using exact 



relativistic scattering and related form factor techniques [|T^ [18] 



In this paper, we are interested in exact off-shell results for two coupled conformal field 
theories^ for which the inter-layer coupling preserves integrability. The on-shell dynamics 
of these models were studied in |20|, 

Let us briefly recall the ideas of |T^, We consider two planar systems corre- 

sponding to two coupled minimal models Aip/p' which interact through a relevant operator 
which leads to an integrable theory. The resulting action can be written : 



A = Mp/pi + Mrj/r,' + A 



p/p' 



12 ^12 



or 



A = Mp/pf + Mr>/v' + A 



p/p' 



■21 ^21 ! 



(2) 



where we denote respectively $12'' ("^21^) ^"^^ "^12^ ("^210 two specific primary operators 
of each unperturbed minimal models and where the parameters A and A characterize 
the strength of the interaction. Here, we will be interested in exact one-point functions 
in such system. 

In section 2 we introduce the notations and those known results which are useful for our 
purpose. In section 3 using the exact result for the VEV of the 6*2^'' ATFT derived in |13[ , 
we deduce the exact VEV < ^l^J {x)^^^l,{x) > for any values of (r, s), (r', s') in the model 
with action (|I|). To do it, we relate the parameter A in (|I|) to the masses of the particles 
and we perform the QG restriction of the C2^^ ATFT with imaginary coupling which leads 



(2)/^{2)> 



to the model For 



'1,2) this VEV can be calculated exactly as well as 



the bulk free energy. The specific case A^3/4 + M.3/A coupled by $12^ $12^ is considered in 
details. It corresponds to two layer Ising models coupled by their magnetization operator 
The previous approach is also extended to the model described by action (^. 
In section 4, we extract some limited information about the asymptotics of two-point 
correlation functions between any pairs of primary operators which belong to the same 



^In literature, the first example of such integrable coupled models was studied in m. 
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or different unperturbed models. More precisely, we will distinguish four different cases 



[a) 
(b) 
(c) 
id) 



< ^l^J{x)^^^l{y) > ior \x-y\ 

< <l>i'^{x)<l>%{y) > for \x-y\ 

< ^^^}{x)^%{y) > for \x-y\ 



0, 



oo, 



and i e {1,2}, 
oo and i G {1, 2} 



which are depicted in figure 1. We finally give some numerical results for various examples 
of coupled minimal models such as two energy-energy coupled tricritical Ising, two coupled 
A5-RSOS models and two energy-energy coupled 3-state Potts model. Perspective and 
conclusions follow in this final section. 




Figure 1 - Two coupled two dimensional models. Short distance 
results are obtained by taking the limit e ^ 0. 



2 Coupled minimal models as restricted €2^^ ATFT 

The ATFT with real coupling b associated with the affine Lie algebra C2^^ is described 
by the action in the Euclidean space : 

A= f d''x\—(d.^? + /i'e-2*^i + fi'e^'^^ + /ie''(^i-^2)l (3) 
J '-8% -I 

where we choosed the convention that the length squared of the long roots is four. As the 
different vertex operators do not renormalize in the same way, we introduced two scale 
parameters /i and n'. The fields in eq. (||) are normalized such that: 

< (^i{x)Lpj{y) >= -5ij In \x - y\^. (4) 

This model possess two fundamental particles with mass Ma which depends on one pa- 
rameter m : 
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TTQj 

Ma = 2msm{—) for a =1,2 (5) 
H 

where we introduced the "deformed" Coxeter^ number H = h{l — B) + B with 
B = and h = = Q are respectively the Coxeter and dual Coxeter numbers. The 
exact relation between m, /i and /i' was found in |]12[ and is given by : 



In ATFT approach to perturbed CFT, one usually identifies the perturbation with 
the affine extension of the Lie algebra Q. Instead, the perturbation will be associated 
here with the standard (length 2) root of C'^^ Removing the last term in the action 
leaves a = 5*0 (4) = SU{2) ® SU{2) model, i.e. two decoupled Liouville models. To 
associate the two first terms of the Toda potential to two decoupled conformal field 
theories, we first introduce for each one a background charge at infinity. Then, the total 
stress-energy tensor T{z) is written T{z) = T^^\z) + T^'^\z) where : 

TW(^) = -^{dip.f + Q,d\., for I G {1, 2} (7) 

ensures the local conformal invariance of the D2 model for the choice Q2 = —Qi = 6+1/26. 
With our conventions 0, the exponential fields 

ga,^,W iG{l,2} (8) 

are spinless conformal primary fields of each Liouville model with conformal dimensions 

A(e">'^»(-)) = -^ + aiQi and i e {1,2}. (9) 

In particular, the exponential fields e"^'"^^ and e^'"^^ have conformal dimensions 1. As 
is well known, the "minimal model" M.-pjpi with central charge c = 1 — 6 can be 

obtained from the Liouville case. Consequently, the D2 CFT can be identified with two 
decoupled minimal models by the substitution : 

h-^i(5, -/i, /i' -/i', (10) 

and the choice : 

/32 = ^2 ^ p/2p' or 0^ = 0i= p'/2p with p < p' . (11) 

Similarly, the primary operators of each minimal model M.p/pi are obtained through the 
substitution aj — > ir]j for j = 1,2 in (^. With these substitutions, the conformal dimen- 
sion of the perturbing operator becomes : 

\,rt = A(e^'^('^^-'^^)) = 3/52 - 1. (12) 



^Differently to ttie simply laced case for which the mass ratios take the classical values, the mass ratios 
for non-simply laced case get quantum corrections. 

•^Here, the length of the longest root is chosen to be 4. 
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As long as we consider a relevant perturbation, we are restricted to choose < 2/3. 
In the following we will consider only the cases for which this condition is satisfied, in 
particular (3"^ = I3\. 

We define respectively {'J'rs''} ^'^d the two sets of primary fields with con- 

formal dimensions : 

A^, = ''^ ^ ~ for 1 < r < p, 1 < s < and p< p' . (13) 

They are simply related to the vertex operators of each minimal model through the 
relation: 

= iV«-i with = -V2 = " (1 - s)f3, (14) 

where we have introduced the normalization factors Nj^^J for each model. These numerical 
factors depend on the normalization of the primary fields. Here, they are chosen in such 
a way that the primary fields satisfy the conformal normalization condition : 

<$«(a;)$«(?/)>cFT= r— Irk;: for ^ E {1,2}. (15) 

For further convenience, we write these coefficients A^^^*-* = N^^\r]^^) where : 

-r(2/32 + 2r//3)r(l/2/32 - r]/P)r{2 - 2(3^)r{2 - 1/2/3^ 



r(2 - 2/32 - 27](3)T{2 - l/2/?2 + r]/(3)T{2(3^)T{l/2(3^) 



iV«(r7) = [-nfi'-f{-2p^) 

and N'^^\7]) = N^^\-r]). 

For imaginary values of the coupling b = the 6*2^^ ATFT possesses complex soliton 
solutions which interpolate between the degenerate vacua. This QFT possesses the QG 
symmetry associated to Uq{Df^) - as we will recall in the next section. In [^] the S- 
matrix of the b'^^ ATFT was constructed using the Uqi^Af^) QG symmetry of this QFT. 
In particular, using the bootstrap procedure, the authors deduced the breather-breather 
5'-matrix. It was also shown that a breather-particle identification holds by comparing the 
5'-matrix elements of the lowest-breathers (breathers with lowest mass) with the S'-matrix 
elements for the quantum particles in real ATFT. Following the conventions of [|l^, the 
particle spectrum in real -82^'' ATFT is given by Mi = 2msin(7r/if) and M2 = m. Using 
the results of [031 we find the identification : 



m = 2Msin(^^) with ^ = (16) 

where we denote M as the mass of the lowest kink. Eq. (|T^) holds for our case due to 
the identification B^^^ = ^ and A^'^ = Df\ 
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3 Expectation values in coupled minimal models 



In ref. [jT2|, |T3| we derived the exact VEV G{a) =< exp{a.cp){x) > for all non-simply 
laced ATFTs using the so-called "reflection relations" which relate different fields with the 
same quantum numbers. We refer the reader to these papers for more details. Although 
the model (|lD is very different from the C^"^ ATFT (^) in its physical content (the model 
([I|) contains solitons and excited solitons), there are good reasons to believe that the 
expectation values obtained in the real coupling case provide also the expectation values 
for imaginary coupling. The calculation of the VEVs in both cases {b real or imaginary) 
within the standard perturbation theory agree through the identification b = i(3. Following 
the analysis done for the Bullough-Dodd model it is then straightforward to obtain 
the VEV of primary operators which belong to different minimal models. With the 
substitutions (|10D one gets Q{r]) = G{ir}) : 



Gil) 



7rV/i'7( 



X exp 



with 



^ dt / 



(i+S) 
2-5 



(A1+A2) 



7r/i7( 



1 + r 



1 + 



t Vsinh((l + ^)t) sinh(2tO sinh((4 - 2^^) 



(17) 



sinh((r7i - r]2)f3{l + ^t) sinh(((r7i - 7^2)^(1 + + 2^ - 2)t) 



sinh(t) cosh(t^) 



+ sinh(2772/5(l + 0^) sinh((2r/2/5(l + - 2)t) 
+ sinh(2r7i/5(l + ^t) smh{{2r]i(3{l + + 2)t) 



sinh((l-Ot). 



As usual we denote 7(0;) = r(x)/r(l — x). The integral in (|r^) is convergent if : 







< (^1 + V2)P < 



;i+o 



and 



1 <{r]i- r]2)(3 < 



3-e 

(1 + 0' 



and is defined via analytic continuation outside this domain. In ( p^Tf ) we defined Ai 
and A2, the conformal dimensions in the imaginary coupling case by Ai = A(e*''^'^^) = 

^ + ^{^^) and similarly for A2 with the change rji — >■ —772- 

If we want to express the final result for the VEV in terms of the parameter A in the 
action (P, we need the exact relation between A and the parameters n, n' in the 6*2^'* 
ATFT with imaginary coupling. We obtain : 



A 



vr/i/i 



(4/?2 _ 1)2 

which corresponds to A = —j2N^^\(3)N^'^\—(3) 



7(4/3V(l-2/5^ 



(19) 
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Like any other ATFT, the model for imaginary values of the coupling has non-local 
conserved charges {Qk, Qk} fo^ k = 0,1,2 generated respectively by the purely chiral and 
anti-chiral components : 



J^y (z) 



and 



(20) 



where the fundamental vector field (p{z,z) = ip{z) +^(^) and {e^}, /c = 0, 1,2 is the set 
of dual simple roots of the non-simply laced affine Lie algebra C2^\ We also define the 
topological charge : 



Hk = J d'^x el.d^ip{x,t). 



Using the equal-time braiding relations for all x,y : 

Jev{x,t)J_e]'{y,t) = q^'"''J_^v{y,t)Jev{x,t) with qi = e I'^fti"" 



2iTr 



(21) 



(22) 



where aki = denotes the extended Cartan matrix of C2^\ one can show that the 

charges {Qk,Qky ^k} for k = 0,1,2 satisfy the quantum universal enveloping algebra 

(2) — 

Uq{D\ ) with deformation parameter q = qq = e ^ . If we express these generators in 
terms of the standard Chevalley basis {E'^, , Hk} by : 



Qk = EU""' and = E^q 



for k = 0,1,2, 



(23) 



then we have : 



[Hk,Hi]=0, 
[Hk,E^] = ±ai^kE^, 



(24) 



k^Ei 



-2H, 



Qi - Qi 



-1 



The Uq{Df^) has two subalgebras Uq{D2) as well as Uq2{C2) where the subalgebra 



Uq{D2) is generated by {Qo, Qq, Q2, Q21 Hq, -^2}- As was found in p3[ the S'-matrix in the 
unrestricted form acts as an intertwinner on the modules of the Uq{D2) representations : 



Pb 



(25) 



where Vp is the module with highest weight p. The scattering of two solitons of species a 
and b is then described by Sa,h{0i2) with 6*12 being their relative rapidity. There are two 



such fundamental multiplets denoted {4} and {6} in ref. In addition, there are also 
scalar bound states and excited solitons depending on the values of {p, p') chosen. 

To understand the restricted Cg^^ (denoted RC'^^ below), we use the general framework 
of superselection sectors (for details see ^). The model (|I|) is a perturbation of 
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the two minimal models by the operator $^2^ $^2^. Each minimal model contains a finite 
number of primary fields (|13],|14D. Using the superselection criterion for the present case0: 

we find j + j € Z where jj denote the representations of Uq{D2) with j the spin-j 
representation of SU{2) with dimension 2j + 1 (and similarly for j). If g is a root of 
unity i.e if eq. (pUD is satisfied then using the Coulomb gas representation condition (|13D 
l<2j + l<p— 1 and l<2j + l<p— 1 one has : 

< j <p/2-l and 0<j<j9/2-l. (27) 
The superselection sectors Ti-^^^ of the RC^^^ model (|l|) are thus : 

T^^^^'' = ^ n^^^^'^ with J + j G Z. (28) 

(i,j)6{0,l/2,...,p/2-l} 



As shown in [20, 21] for the unitary series {p' = p + 1), after the quantum group retriction 
( [Tl| ) for p > 3 the fundamental solitons in the {4} and {6} representation of Uq{D2) 
become the RSOS kink K-'.-'-. . - . These kinks interpolate between different vacua |0,- > 

3232 3\3\ •'^■'^ 

and 10^2^2 > which are connected using the Uqi^sh) fusion ring at g = — e~*p: 

mm(ji+j,p-2-ji-j) 

Jixj= ^2 (29) 

32=\3i-3\ 

and similarly for j. However, for p = 3, (j, j) G {0, 1/2} and then the {6} is projected out 
of the spectrum, leaving only the {4}. We refer the reader to ^ for more details. 

From the previous remarks and the identification D2 = SO (4) = SU{2) (g) SU{2), 
by analogy with ^ the primary fields ^^i] and ^f) commute with the generators in 
(H) (for k = 0,2) of the subalgebra Uq{D2) C Uqiof^). If one interpret the fields 



$2j+i i(-2)*^2j+i i'-'^'' highest component fields in the multiplet, it can be shown that 

the primary operators and are not invariant with respect to Uq{D2). Together 
with some non-local fields they form finite-dimensional representation of this algebra. 
Consequently, the VEV in (|T]) should take into account the factor d^^ coming from 
the QG restriction of the QFT (|^). Following the conjecture of it takes the form : 

sm(^^^2l±ll\p'r -ps\) sin(^i^bV - ps'l) 
sinC-^ip'-p)) sin(^(p'-p)) • ^ ^ 



*We only consider the holomorphic part of the primary operators but keep the same notation. 
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Using the notations introduced in the previous section, and eqs. (0), (IT^); outcome 
for the VEV between different primary operators is : 



4-2g 



< 0,;.|<l>lV(x)<l>l'L(x)|0,:. > 



IV rs 



Mi 



-7rA7(T^)(l + 0~ 



(1+0 
2-€ 



(A,,+A,,,^,) 



X exp Qi2((l + 0^ - 2^5, (1 + 0^' - 2^5')- (31) 
The function Qi2(^, d') for |0 ± < 4^ and ^ > is given by the integral : 



Qi2ie,e') = 

with 



°^ dt / 



t Vsinh((l + i)t) sinh(2tO sinh((4 - 2^^) 



02 + ^'2 -2(1-^)2 



-2t 



cosh{{e + e')t) cosh{{e - e')t) 



cosh((2 - 2Cjt) 



sinh((l-Ot) cosh((4-20t) 



+ 



cosh((0 + e')t) + cosh((0 - e')t) - cosh((2 - 2^^) - 1 



sinh(t) cosh(t,^). 



and defined by analytic continuation outside this domain. Notice that eq. (jST]) satisfies: 



(1) 



.(2) 



> 



(32) 



A particular case of eq. (|3l|) is the expectation value of the perturbing operator which 
can be calculated explicitly to give the result : 



< 0-|$i2^(x)$f2Hx)|0- > 



33 



-7rA7(Yq:|) 



7(fS)7(M) 



^2^(1 + 07(4^^)7(4^) 



<i - 2) 



7( 



4-25' 



.(33) 



By using the exact relation between the mass parameter m and the mass of the kink M 
(1T6| ) and eqs. (||), ([T9| ) we immediately obtain the relation between M and A : 



M 



4-25^ V4_25' 



4-25 ' 



4-2{ 



(34) 



Consequently, according to eqs. (|l^), (|33|) , (^) and /52 < 2/3, the perturbed CFTs (|l]) 
develop a massive spectrum for (Jm(A) = 0) : 



0<^<l/3, A>0 i.e. < - < -, 

p' 2 



(35) 



(u) 1/3 < ^ < 1, 



1 P 

A < i.e. - < ^ < 1. 
2 



where ^ = -p ■ particular, the condition {ii) is always satisfied for the coupled unitary 
minimal models defined by (|l]). 
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Finally, the expectation value (^) can be used to derive the bulk free energy : 



12 



— lim — In Z, 

V^ooV 



(36) 



where V is the volume of the 2D space and Z is the singular part of the partition function 
associated with action (|T]). Using 



dxf: 



12 



and eqs. 



the result for the bulk free energy follows 



sin( 



12 



4-25 



) sin(^) 



4-2^' 



sm 



■ vr(l+g) ' 
. 4-2g ' 



(37) 



(38) 



As was suggested in a particular case of the model (|l|) can be related to the 

ATFT with imaginary coupling. 



SG theory at the reflectionless point as well as the Dg 
Let us first consider the SG theory with action : 



(1) 



A, 



SG 



16n 



(39) 



This model is integrable and its on-shell solutions, i.e. spectrum of particles and S'-matrix 
are well-kown . This theory contains soliton, antisoliton and soliton-antisoliton bound- 



state (breathers) -B„,n=l,...,l/^ where ^ = -r^j^- The lightest bound-state Bi coincides 
with the field in the perturbative treatment of the QFT (pOj). Its mass is given by 
mi = 2M5G'sin(^) where Msg denotes the soliton mass of the SG model. For p = 3, 
p' = 4 the model (||) corresponds to two magnetically coupled Ising models. It is a c = 1 
CFT perturbed by the operator $12^ $12^ with conformal dimension 1/8. Indeed, the VEV 
of this operator must be simply related to the VEV of < cos(/30) > in the SG model for 
= 1/8, which is also a c = 1 CFT with perturbing operator of conformal dimension 
1/8. For this latter value, ^ = 1/7 and the SG model possesses 6 neutral excitations : 



a = 2MsG sin 



14 



a = 1, 6 ; 
mi < 1712 < Msg < ^3 < ■■■ < "^6- 



(40) 



Similarly, for = 3/8, ^ = 3/5 the spectrum of the model (|I]) was described in and 
summarize here. Here, we denote the lowest mass by M. If we compare the spectrum of 



each model, using the notations of I^S] we have : 

mi = M, m2=m (1), MsG = 'm (2), etc... 



Furthermore, the bulk free energy of the SG model (|H^) is well known 
by : 



SG 



(41) 

and given 

(42) 
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By expressing it in terms of the lightest particle mi using (^) and the relation 

8 sin(7r/14) sin(37r/14) cos(7r/7) = 1 we find the same result as in eq. (|38| ) for = 3/^ 

i.e. ^ = 3/5. Also, using the results of for = 1/8, one has : 



A 



27(1/^ 



MsGv^ r(4/7) 



TT L 2 r(l/14) 
Similarly, the VEV of the perturbing operator is given by 



(43) 



,^ _ 8 r(i/7)r(6/7) 

< cos(/^0j >42=i/8- 7r2(4/7)r2(i3/i4) 



TT 



7(1/^ 



8/7 



2^ ■ 



(44) 



Comparing eqs. (|43|), (|4^ ) with eqs. ( |33D and (|^) respectively for ^ = 3/5 and using eq. 
(|4l|), one find the relations : 



A = 22 A and 



COs(/3(/))|^2^1/g 



(45) 



in perfect agreement with eq. (B.21) in 



Let us now turn to the relation with the D^^ ATFT. The calculation of the bulk free 



energy for the simply-laced Dg ATFT gives p8| 

2 

rr 



-^8 lb 



8 sin(^)sin( 



(46) 



where B is defined in section 2 and the mass of the particles are related with the mass 



2m^, ml 



sm"l ^1 for a 



parameter m by : mg = rrij 
coupling, one expects the particle-breather identification : 



2—2 ^i„2fna 
14 



8m^ sin^l 



14' 



2Msin(g))' 



1, 6. For imaginary 
(47) 



where ^ is defined as in eq. (|16D. This yields : 

M2 



sm 



, 14. 



16 sin(fL)sin(^) 



(4J 



which for ^ = 7 is in perfect agreement with eq. ( P^ evaluated at ^ = 3/5. 

It is also interesting to study the behaviour of the model (|l|) for p' = p + 1 in the limit 
p — >■ cxo. If we consider the sine-Gordon model with action ( pUD there is an equivalent 
description in terms of the massive Thirring model |3^ : 



A 



MT 



(49) 



At (3'^ = 1/8, the operator can be expressed [|| [29) in terms of two independent Ising spin fields, 

12 



cr(») = for i^l,2 
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where ip,ip is a Dirac field. The four fermion coupling constant g relates to P in (^) by 
^ = ^ — 1. Also the free fermion point is reached for — |. Using the results of ref. 

Q concerning the one-point function < e*"*^ > in the SG model one gets : 

Mqr 1 

< cos(/30) >^ lim MsGT{e) and A ^ — for -. (50) 

e->0 TT 2 

Using the boson-fermion correspondence, one has the identification ipip = cos(/30)/7r and 
< ipip >— lim ^^^^^r(e). Let us now return to the case p' = p + 1 and the limit p ^ oo 

in (|l|). It corresponds to — ^ ^, i.e. ^ — 1 and similarly to eq. (|50|) the VEV of 

the perturbing operator $12^ $12^ also becomes singular. Consequently, the perturbing 
operator in (|1]) can be conveniently rewritten in terms of two Dirac fields (with the same 
mass Mf) of two independent massive-Thirring models at the free fermion point : 

A$S2^$Sf ^ -Mf^tjj - Mf^'ijj'. (51) 

In fact, this situation is not really surprising : when p ^ 00 the coset algebra of each 
Mp/p+i reduces to a level-1 SU{2) current algebra. The two models are coupled via their 
primary fields in the spin | representation of dimension ^ [^. As expected, in this limit 



the model becomes the free field theory of two complex massive fermions with undeformed 
SO {A) symmetry. 

Finally, let us consider the case p' = p + 1 and p = 2 in the action (|l]). Each minimal 
model is then identified to A^2/3 with central charge c = 0. For each model, the unitary 
representation is the vacuum one with conformal weight An = A12 = 0, for which cases 
we have the identifications $[2 = I*^*^ fo^ * ^ {!> 2}. This corresponds to the choice = ^, 
i.e. ^ = I and, using eqs. (0), (0), (^1), one can check that 

< 0|<I>[2^<I>S2^|0 >= 1 and /12 = A, (52) 

as expected. 



Let us now turn to the model associated with action (|^). Due to eq. ([Tl| ) it corresponds 
to the choice = The condition ip > 3p' guarantees that the perturbing operator is 
relevant. Then, the vacuum structure is expected to be similar to that of (|l]). From eqs. 
and ( PID and the substitutions : 



p ^ p', (r, r') • 
the result for the VEV immediatly follows : 



s, s 



1 + e 
3^-1 



< 0,5|<l>lt)(x)<l>l?!,(x)|0,5 > 



-7rA7(^)(20 — 



7(|)7(^) 
X exp Q2i((l + Or - 2^s, (1 + Or' - 2Cs') 



(53) 



(54) 
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The function Q2i{0,0') is given by the integral : 
dtf "^21(0, 6', t) 



with 



— - 



t Vsinh((l - ^)t) sinh(2tO sinh((3 - d^t) 



-2t 



^2i(^,^',t) = cosh((^ + ^')^)cosh((0-^')^) -cosh((2-20t) 
X sinh((l - ^)t) cosh((3 - 5^)t) 

cosh((^ + 9')t) + cosh((0 - 9')t) - cosh((2 - 2^^) - 1 
X sinh((3^ - l)t/2) cosh((l + 0^/2) 

and defined by analytic continuation outside this domain. The prefactor associated with 
the QG restriction dHy^, = dll ,.,^,\p^p/ . For (r, s) = (r', s') = (2, 1) becomes : 



<0^..|$«(x)$g)(x)|0^.-.> 



7(|)7( 



1+5 



3$-l 



TT 



(3 - 50 7(#3) 



(55) 



with the relation between the mass of the lightest kink M and A : 



4? 



7(i)7(^) 



2€ 



55-3 



(56) 



For the coupled minimal models defined by (0), the massive phase corresponds to the 
domain : 



{m) - < ^ < 1, 
5 



Q 

A<0 i.e. - < ^ < 1. 

4 p' 



(57) 



One also obtains the bulk free energy associated with action 

f21- 



^^^M^)sin(^: 



sm 



(5J 



4 Application and concluding remarks 



Accepting the conjectures (pT]) and (p^D, one can easily deduce interesting predictions^ 
for the short and long distance asymptotic of two-point correlation functions in the model 
([1|) or (0). For each case depicted in figure 1, we can express the result in terms of : 



<0^.-.|<l>«(x)<l>^?i,(x)|0,.>=C,,,^^ 



(59) 



^The same method was applied in Q to obtain a prediction about the long distance asymtotic of 
two-point correlation function < (Tocr„ >„^oo in the XXZ spin chain. 
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given by eq. ( PT] ) or (0), respectively. First, using the short distance approximation 
(operator product expansion) and eq.(|) we get e^'^'^i^^^e^^'^^Cy) ^ ^i(aM^)+bMy)) _ Then, 
using the Coulomb gas representation of each primary operator (|14D one has 

< 0^.-.|<l>«(x)$S,(y)|0^.-. > — < 0^.-.|$«(a;)<l>S,(a;)|0^.-. > . (60) 

\x—y\^0 

In this limit - case (a) - , the two-point functions then become : 

<0f^\<l>g\x)<l>%{y)\0f^> ^ gg,,,,. (61) 

Secondly, in the long distance approximation the asymptotic two-point function simply 
reduces to the product of two one-point functions as < e^'^^'^'^^^e*'' '^'^^^ >— e*'' '^*^^'^ >< 
Q^b.ifiiy) ^ -vv^Jien \x — y\ ^ oo. Then we obtain - case (6) and (d) - 



< Of,\<^^r's\^)0Ay)\Of, > < o,-.|$«(x)i(2)|o.. X o.|i(^)$S(x)|o.. >; 

\x—y\ — >oo 

< 0,j|$«(x)$S(y)|0^.-. > ^^-^^ < 0^.-.|<|.«(x)l(^)|0^.-. X 0.|<|.S(x)l(2)|0^.-. > . 
Indeed, it gives : 

<0^-^\^i'J{x)^%{y)\0^.> — Giu\Sng^s',n + S^2gir's] for {1,2} (62) 

with 6iir the Kronecker symbol. Obsviously similar results are obtained for the two-point 
function < 0--.m'}{x)^l%{y)\0 f- > using the Z2 symmetry (1 2). 

Finally, using the fusion rules in the short distance approximation (the two pri- 
mary fields belong to the same space of states) x <l>^!^^, — > $^!/^;/, the two-point 
function is expanded in terms of the one-point functions < Ojj\^^j^}g„{x)I^'^^Ojj > or 
< Qjj\I^^'^^^^}g„{x)\Ojj > for i G {1,2} respectively. Then, we have for i = 1 - case 

< OfM'^ix)0Ay)\Of, > , — > E - (63) 

and similarly for i = 2 where Cl!Ar's' ^^e the structure constants of the minimal model 



operator algebra |33 



Example 1 : Two magnetically coupled Ising models 

It is now straightforward to compute different two-point correlation functions in one of the 
simplest (non-trivial) cases : two-magnetically coupled Ising models in the massive phase. 
It corresponds to = 3/8 i.e. ^ = 3/5 in (|31|). In this case we have the identification 
^12 = A"^"^ with conformal dimension A^- = 1/16 - the spin operator - and $^3 = e^'-* with 
conformal dimension = 1/2 - the energy operator. There are two degenerate ground 
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states |2y, |2T| denoted |Ooo >= |— > and |0i/2 1/2 >= |+ >• For simplicity, we write 
< ±|...|± >=< ... >±. Sometimes, the reflection relation : 

G{m,V2) = Sl{v2)G{vi, -V2 + l/P) (64) 

with 



Sl{V2)= - 7r/i'7(-2/32 



X 



r(2 + 2ri2p - 2/32)r(2 - r]2/P - 1/2(3^ 



is useful for analytic continuation. Using eq. ( ^I]) and (in, is = df^ n for all s we obtain 
for instance : 

< a« >^ = ^±^^2 = ±1.297197220. ..(-A)^/^^ 

< e« >^ = ^±^^3 = -2.278284275. ..(-A)^/^; 

< (T(^Ho)fT^^nO) >± = ^i1,i2 = 1.698928047... (-A)^/^; 

< aW(0)(T(')(oo) >± = {Gti,i2f = 1.682720628... (-A)^/^; 

< a(^)(0)e('H0) >± = e?i^2,i3 = ±3.311880669.. .(-A)^/^^ 

< (T(^)(0)e(2)(oo) >± = g^i^i 12^^1^1,13 = ±2.955384028... (-A)^/^^ 

< e(^)(0)e(')(oo) >± = 13)2 = 5.160349412.. .(-A)^/'; 

where the parameter A is related to the mass of the lowest kink by : 

A = -0.2379062104.. .M^/^ (65) 

Notice that < a^^\Q)a^'^\Q) >± and < a^^\Q)a^^\oo) >± differ by less than 0.7% as 
expected []. 

Example 2 : Two energy-energy coupled tricritical Ising models 

The case p = 4, p' = 5 in (|1|) describes two tricritical Ising models which interact 
through their leading energy density operators $^2 = ^ of conformal dimension = 
1/10. It corresponds to = 2/5 i.e. ^ = 2/3 in (|3l|) . Beside e^^^ and the identity 
operator, each minimal model contains the sub-leading energy density operator = e'*-*^ 
with Ae' = 3/5 ("vacancy operator") , two magnetic operators $22 = '^''^^ with A^- = 3/80, 
$21 = cr'*''^ with Ao-' = 7/16 and $^4. Due to the obvious property d^^ir's' — '^11 T-'s''^rlii) 

^For instance, it is known that form factors are able to reproduce with high accuracy the UV behaviour 
of the correlation functions IH, pi], p2| . 



15 



similarly to the previous case we find for instance for any vacuum \jj > : 
<a«>. = ^g^,, = 4,, xl.l44656674...(-A)='/^^ 
< e(i) > = ^^2,11 = <2,ii X 1.529866659... (-A)^/^ 
<a«(0)a(2)(0)>^.. = 6;i,22 = 4,22xl.315726811...(-A)3/32; 
< a(i)(0)a(2)(oo) >^.. = gg,^i,2 = 4.22 X 1.310238901...(-A)3/32; 

<6«(0)e(2)(0)>^.3 = ^|i2 = 4i2 x 2.419476973...(-A)^/^; 
< e(^)(0)e(2)(oo) >^.^. = gg,^i,2 = 4,12 X 2.340491994...(-A)i/^ 
where the parameter A is related to the mass of the lowest kink by : 

A = -0.2566343706.. .M^/^ (66) 

Notice that < (T(i)(0)a(2)(0) >■■ and < a^^\0)a^^\oo) < e^^\0)e^^\0) and < 
e*^^^(0)e*^^^(oo) >jj differ by less than 0.5% and 3% respectively. 

Example 3 : Two coupled RSOS models 

The case p = 5, = 6 in (|T|) describes two RSOS models coupled by their primary 
operators with conformal dimension A12 = 1/8. It corresponds to = 5/12 i.e. 
^ = 5/7 in (pTI). Each minimal model also contains the primary operator $22 with 



A22 = 1/40. As before we obtain : 

< *S >j] = = 4,11 X 1.090446894... (-A)i/=^0; 

< *S >j] = <3iiii = 4,11 X 1.726352342... (-A)^/^; 
<$g(0)<l>g(0)>^.^. = 6^^2,22 = 4,22 X 1.191588988... (-A)i/i^ 

< $g(0)$g(oo) = gi,iigil,22 = 4,22 X 1.189074429... (-A)i/i^; 
where the parameter A is related to the mass of the lowest kink by : 

A = -0.2697511940.. .M3/2_ (g7) 

Notice that < $S(0)$22 (0) > - and < $^2^(0) $^2^ (00) > - differ by ~ 2%. 

Example 4 : Two energy-energy coupled 3-state Potts models 

The case p = 5, = 6 in (||) describes two 3-state Potts models coupled [l^ by their 
energy density operator $21 = with conformal dimension A21 = 2/5. It corresponds 
to ^ = 5/7 in (54). Each minimal model also contains the primary operator $23 = " 



the spin operator - with A23 = 1/15. We obtain for instance : 

<a«>. = ^tn = 4,iiXl.9079...(-A)i/3; 
<a« (0)^(2) (0)>,^ = ^i,23 = 4,23 x 4.50...(-A)2/^ 
<a^'\0)a^'\oo)>^. = ^4n^a23 = 4,23 x 3.64...(-A)2/3; 
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where the parameter A is related to the mass of the lowest kink by 



A = -0.2612863655.. .M^/^ 



(6J 



To conclude, we would like to mention that here we studied only a special case of a 
much more general class of integrable coupled models [|19], ^ There exist many other 
examples which can be worked out along the same lines. Let us also note that the exact 
form factor techniques as well as the truncated conformal space approach may be used 
and similar numerical analyses can be performed for the correlation functions. 
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